THE UNIVERISTV OF ALABAMA IN HUNTSVILLE 


Analysis of SOFCAL Calibration Data 


FINAL REPORT 


Submitted by: 


Dr. John C. Gregory 
Department of Chemistry 
The University of Alabama in Huntsville 
Huntsville, AL 35899 
(205) 895-6028 


December, 1993 


(NASA-Ca-193910) ANALYSIS OF N94-25000 

SOFCAL CALIBRATION D ATA Final 
Report (Alabama Univ.) 20 p 

Unci as 


G 3/35 0206777 


INTRODUCTION 


Supporting Studies for A Scintillating Optical Fiber Ionization Calorimeter (SOFIC) 

/ 

Ionization calorimeters determine the total energy of a particle by absorbing within 
the mass of the calorimeter the entire energy of the particle or a significant (and 
determinable) fraction of it. The JACEE group has used passive ionization calorimetery, 
employing photographic emulsions as the detector medium, to measure the charge 
composition and energy spectrum of cosmic rays up to, and exceeding, 10 14 eV. The 
SOFIC approach depends similarly on the use of three-dimensional nuclear-electromagnetic 
shower theory to relate the ionization deposits obtained during a shower to the energy of 
the particle causing the shower. Bundles of thousands of scintillating optical fibers are read 
out using image-intensified CCD’s. Such an event may be triggered by a fast shower 
detector placed under the instrument or if only heavy particles are of interest, from a fast 
primary Cerention detector placed above the calorimeter. 

In the first study, fast Hammatsu photomultiplier tubes were purchased and tested 
for possible application for a triggering purpose. 

In the second study, some refinements have been made to the theoretical treatment 
of hadronic interactions in the central collision region. These will be helpful in improving 
the simulations necessary for observations of high energy cosmic ray nuclei with a 
SOFIC. 
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1.0 Introduction 


The Hamamatsu RA1250 is a large (5 inch photo-cathode) photo-multiplier, which 
employs a liner focused 14 stage dynode geometry to produce fast (t_rise - 2.2 ns) 
response. The entrance window is borosilicate giving a spectral response from 300 to 
650 nm. The quoted quantum efficiency of the Bialkali photo cathode is 22%. The 
results of measurements made on ten of these tubes with novel base electronic is 
presented. 


2.0 Novel Base electronics 

Photo-multipliers are conventional operated with a passive resistor base supplying 
appropriate bias voltage. Some experiments have been made with electronic bases which 
can offer low power consumption (ref 1), or high count rate capability (ref 2). In these 
measurements we have experimented with an active base using FET. The base circuit 
(ref 3, Fig 1), shows how the dynode chain is biased with FET, which also allow a 
slowly varying signal to be extracted via the dynode path as well as a fast anode 
signal. Such an arrangement is ideal suited to applications where fast, and slowly varying 
signals need to be measured with the same photo-multiplier. The dynode circuit performed 
satisfactorily. 

3.0 Measured Performance 


3.1 Linearity 

The HM photo-tubes were stimulated with a light pulsar (241Am in Nal). The response as 
a function of voltage was linear as shown in Fig 2. 

3.2 Relative gains 

Although the tubes were purchased as a matched set there was a gain difference of two 
between the best and worst tubes (Fig 3). 

3.3 Position response 

The signal variation with photo cathode position was investigated by recording the photo 
tube output with the light pulsar placed at different positions across its surface (Fig 


3.4 After pulsing 

Most of the tubes displayed after pulsing. Fig 5, and Fig 6 show examples of the best and 
worst Five of the tubes were exchanged with Hamamatsu, but the replacements all had 
some after pulsing. It appeared that there is a problem which Hamamatsu could not 
correct. 


3.5 Photo-tube oscillation 


The Hamamatsu tubes were placed in a Lambertian sphere, and stimulated with a UV 
laser. Wavelength shifting aerials were employed to wavelength shift the light in to the 
tubes acceptance window. Typical photo tubes response is shown in Fig 7 The 
oscillatory structure shown is inherent to the photo tubes, and not associated with the 
FET base electronics: oscillations were still visible when the FET base was replaced 
with a passive resistor base constructed according to the prescription provided by 
Hamamatsu. Similar structure was also observed, with no external stimulus, when the 
tube was occasionally struck with a cosmic ray muon. It appears this tube design has a 
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Figure captions: 

Fig 1. The FET base bias circuit 

Fig 2. The output response of the Hamamatsu RA 1250 photo-tubes as a function of 
bias voltage. The right hand scale refers to the Hamamatsu serial number. 

Fig 3. Relative gain variations of the Hamamatsu RA 1250 photo tubes. 

Fig 4. Position response of the Hamamatsu tube: serial number 37. 

Fig 5. Best after pulsing response; tube number 80 

Fig 6. Worst after pulsing response; tube number 74. 

Fig 7. Typical oscillatory structure observed in the anode output of the RA 1250 tubes 
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Gluon Clusters Formation and Their Subsequent Decay 
Into Hadrons in a High Energy Nuclear Collision 
For a Scintillation-Fiber-Calorimeter 
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Abstract 

When two heavy nuclei collide with extremely high energy, a large number of 
hadrons will be created. In this report, we use a simple thermal model to predict the 
number of gluons produced in the central collision region in a high energy nuclear 
collision. We then use a simple symmetry model to calculate the formation of two-gluon 
and three-gluon clusters and their subsequent decay into mesons and baryon-antibaryon 
pairs. The developed scheme of hadronic interactions will help improving a simulation 
program for observations of cosmic nuclei with a Scintillation-Fiber-Calorimeter. 


I. Introduction 

When two heavy nuclei collide with extremely high energy, it is expected that a part 
of their energy will be earned away by the leading particles which have the baryon number 
of the colliding nuclei, and a large portion of their kinetic energy will be converted into heat 
so that a quark-gluon plasma will be created in the collision region. 1 This plasma 
subsequently hadronizes into baryon-antibaryon pairs and mesons as it cools. This so 
called central rapidity region has been the interest of many recent theoretical and 
experimental studies. 2" 6 This report oudines our recent theoretical study. We use a very 
simple thermal model for the gluon production and a simple cluster model for the formation 
of gluon clusters and their subsequent decay into hadrons. With a simple assumption of 
thermal equilibrium, the number of gluons produced and their energy distribution are 
calculated in terms of the incident energy of the colliding nuclei and their interaction 
volume. The gluons produced subsequently form colorless clusters. The number of two- 
gluon and three-gluon clusters and their mass distribution are calculated from a simple 
symmetry model. The decay of these clusters into baryon-antibaryon pairs and mesons are 
then calculated using an SU(3) flavor symmetry Lagrangian. The particle composition in 
the central rapidity region can be completely predicted without any adjustable parameter and 
can be easily compared with experimental results. This work is still in progress and shall 
be completed very soon. 

II. Gluons Production 

When two relativistic heavy nuclei collide head-on, viewing in the center of mass 
system, we have two highly Lorentz contracted pan cakes passing each other. 1 Shortly 
just before and after they pass through each other in a time of the order of one fermi/c, the 
quarks and gluons in the nuclei are excited and a large portion of their kinetic energy is 
converted into heat within a small cylindrical volume bounded by the two receding nuclei. 
Within the volume, it is expected that a thermal equilibrium is reached and this volume can 
be treated as a black body cavity. Since the gluons are produced more easily than the 
quarks, and they are spin 1 massless particles, they, like the photons, satisfy the photon 
statistics and their energy distribution is governed by the Planck formula 
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U(E,T) = 


( 1 ) 



where 6 — 1/k T. A factor of 8 is included in Eq(l) due to the fact that there are eight kind 
of gluons instead of one kind of photon. Integrating Eq (1) gives the energy density: 


U(T) = 


Sn 2 1 

IT/F 


( 2 ) 


Eq (2) can be used to estimate the temperature of the cavity. Let us assume a 160 GeV/n 
lead beam on lead in the laboratory system, or a total center of mass energy W=3600 GeV. 
Let us assume a third of this energy is carried away by the products of each of the two 
receding nuclei and the remaining one third is left in the central collision region. Therefore, 

UV= 1200 GeV (3) 


where V is the cavity volume. Let the radius of the cavity be the lead nuclear radius, 7 fm„ 
and the length be 2 fm. From Eqs (3) and (2), the temperature of the gluon plasma is 
estimated to be 270 MeV. With the temperature worked out, Eq (1) can be used to calculate 
the gluons average energy and their energy dispersion: 


£ _ fEU(E,T)dE 
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1.03 GeV 


(4) 


and similarly, AEg = 0.55 GeV. 

Eq (4) and Eq (3) tell us that for a 160 GeV/n lead on lead head-on collision, one 
expects that a total of about 1 166 gluons are produced in the central collision region with an 
average energy Eg = 1.03 GeV and AEg = 0.55 GeV. 

III. Clusters Formation 

In the last section, we estimated the number of gluons and their average energy 
produced in a high energy nuclear collision. In this section we study the formation of 
gluon clusters and their properties. Since gluons carry color, they cannot decay into 
hadrons directly, but they may form colorless clusters and then decay into hadrons. We 
expect the most abundant clusters will be formed by low number of gluons, namely, two- 
gluon clusters, Ce and three-gluon clusters, CcoP^ 

Ce : Ce is formed when one gluon is combined with another gluon of opposite 
color, for example a b r gluon is combined with a r b gluon. We may calculate the 
probability of such a formation as follows. Let f be the frequency of one gluon in the 



vicinity of another gluon. Then there is a 1/8 chance it will meet the right gluon and form a 
Ce, and a 7/8 chance it will meet the wrong kind and stay as a single. If it does not meet 
the right gluon on the first time, it may meet the right gluon on the second, or on the third 
time ... So the total probability of forming a Ce is: 


/><c £ ) = I/ + Z/xi /+ (Z / ) xi /+ ... =i X 
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Cco: Cco is formed when one gluon, say br, is combined with two gluons of the 
right color, rg and gb. This calculation is more complicated than the previous one. We 
may proceed as follows. Say a br gluon, if it meets either a rg or a gb on the first 
encounter, it needs to add a gb or a rg on the second encounter or on the third encounter 

.... Hence, the probability of forming a Co) if it has met one of the two right gluons on its 
first encounter is : 


P = -f | 
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But there is also a — chance it meets neither a rg nor a gb on the first time. Out of this - 

1 - 8 
chance, there is a - chance it will meet a rb gluon and leave as a Ce cluster, a - chance 
8 8 


it will meet the others and will stay as a single. If it stays as a single, it always has a 
chance to form a Cco by meeting the right gluons. Hence the total probability of forming a 
Cco is 


P(Cw) = — ^ 
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If we assume that eventually all gluons form colorless clusters and decay into hadrons and 
they only form two-gluon and three-gluon clusters, we have 


P(Ce) + P(C<y) = 1 (7) 

Eq (7) determines / uniquely. Substitute the value of / into Eqs (5) and (6), we find 


P( Ce) = 0.64 and P( Cco) = 0.36 



With 1166 gluons on hand, they will form 316 two-gluon clusters Ce and 178 three-gluon 
clusters Cco. 

The mass distribution of Ce and Cco : Ce is formed by two gluons each having an 
energy Eg and width A Eg. The combined mass of Ce is i/2£« 2 (i -cos#) where 6 is 

the angle between the momentum vectors of the two gluons. Hence the average mass of 
Ce is 


M(Ce) = (^2Eg 2 (l-cos0)} = = 1.37 GeV 

The total width of Ce is partly due to the energy uncertainty of the gluon, AEg, and partly 
due to various possible values of the angle 9 . The part of AM due to the uncertainty of the 

angle 9 can be worked out just like that for M( Ce) and is The total width of Ce 

is therefore: 


A M{Ce) = 


1 ||Eg 2 +2(AEg) 2 =0.92 GeV 


Cco is formed by three gluons. Its average mass and width can be calculated in a 
similar way. The calculation is straight forward but lengthy, one needs to take an average 
over the angles between the three momentum vectors of the gluons. The final results are 

M(Ctu) = 2.47 GeV 
AM(Cto) = 1.06 GeV 

We note here that the widths of Ce and Cco are comparable to their masses and are about 
IGeV which are large enough to prevent them from being observed as resonances and 
small enough for them to be considered as clusters. 

IV. Hadronization 

In this section we study the decay of Ce and Cco into mesons and baryon- 
antibaryon pairs. It is commonly proposed that if gluons combine to form clusters, they 
are likely to be formed with low orbital angular momentum between gluons. This means 
that Ce and Cco shall have the quantum number JPC = O' 1 " 1 " and 1“, respectively. 
Furthermore, they are colorless, quark neutral and flavor singlet. Ce, being made of two 
gluons, can only decay into two mesons. Cco , being made of three gluons, can decay 
either into three mesons or a baryon-antibaryon pair. If we assume an SU(3) flavor 


symmetry, the effective Lagrangrian for various gluon cluster decay processes can be 
written out, for example, for the Ce — > PP and Ce — > W decay processes. 

Lef f = g[ Ce Tr (PP) + CeTr (Vp Vp)] 

where g is the coupling constant, P and Vp are the well known 3x3 SU(3) matrices for 
the pseudo-scalar meson octet and vector meson nonet. Once we have these Lagrangians, 
the branch ratios for various decay processes can be easily worked out. They depend on 
the product of two factors, the SU(3) symmetry factor and the available phase space. 
Knowing the masses of the gluon clusters and the number of these clusters, we can finally 
calculate the hadron composition in the central rapidity region in a high energy nuclear 
collision. 

The last stage of this calculation is still in progress. It shall be completed soon. 

In conclusion, we used a simple thermal model to predict the number of gluons 
produced in the central collision region in a high energy nuclear head-on collision. We then 
used a simple symmetry model to calculate the formation of gluon clusters and their 
subsequent decay into hadrons. 
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